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BAND FUNCTIONS IN THE PRESENCE OF MAGNETIC STEPS 
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Abstract. We complete the analysis of the band functions for two-dimensional 
magnetic Schrodinger operators with piecewise constant magnetic fields. The 
discontinuity of the magnetic field can create edge currents that fiow along the 
discontinuity that have been described by physicists. Properties of these edge 
currents are directly related to the behavior of the band functions. The effec¬ 
tive potential of the fiber operator is an asymmetric double well (eventually 
degenerated) and the analysis of the splitting of the bands incorporates the 
asymmetry. If the magnetic field vanishes, the reduced operator has essential 
spectrum and we provide an explicit description of the band functions located 
below the essential spectrum. For non degenerate magnetic steps, we provide 
an asymptotic expansion of the band functions at infinity. We prove that when 
the ratio of the two magnetic fields is rational, a splitting of the band functions 
occurs and has a natural order, predicted by numerical computations. 
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1. Introduction and motivation 

In this paper, we investigate the spectral properties of magnetic Hamiltonians 
£fa = (-zV - Af = Dl + {Dy - ab{x) f, (1.1) 

in L^(M^) with the following discontinuous magnetic field 

B{x,y) = &illB_(a;) + &2lR+(a::), 

where b = (61,62) S and '^x{x) is the characteristic function on A C K. An 
associated vector potential is given by : 

A{x,y) = (0,a[,(x)), ab(a:) = 6ia;lR_(a;) + b 2 xtw+{.x). 

These operators naturally appear in several models involved in nanophysics such 
as the Ginzburg-Landau model of superconductivity (see [31]) and in quantum 
transport in a bidimensional electron gas. In particular in an unbounded system, 
inhomogenous magnetic fields may play the role of a quantum waveguide ([23, 30]). 
More recently it appears that such magnetic barriers also induce transport and act 
as waveguides in graphene by creating snake states along magnetic discontinuities 
(see [22] and [11] where piecewise constant magnetic fields are considered). These 
physical properties can be described with the corresponding magnetic Schrodinger 
operators that exhibit a variety of interesting spectral and transport properties. In 
particular, there may exist edge currents that propagate along the edge at a: = 0. 


PDH was partially supported by NSF through grant DMS-1103104. 
Version of October 31, 2016. 



2 


P. D. HISLOP, N. POPOFF, N. RAYMOND, AND M. P. SUNDQVIST 


The conductivity corresponding to these edge currents is, in fact, quantized. Most 
of these properties can be deduced from the study of the band functions associ¬ 
ated with the Hamiltonian (see for instance [ 21 , 8]). We will describe later some 
consequences of our analysis of band functions. 

The symmetries of the Hamiltonian allow us to classify various configurations 
of the magnetic field (t>i,62). We notice that = £_(, so that we may assume 
that 62 > 0 . If S' denotes the symmetry {x,y) 1—>■ {—x,—y), we have = 

C-b2,-bi = ^b2,bi- For H > 0 , we introduce the L^-unitary transform 

UBijix,y) = 

and we have 

Ub^C^Ub = B-^CBb- 
These considerations allow the following reductions: 

( 1 ) If bi or 62 is zero we may assume that 61 = 0 and, if 62 ^ 0 , we may assume 
that 62 = 1 - We call the case (61,62) = ( 0 , 1 ) the “magnetic wall”. 

( 2 ) If 61 and 62 have opposite signs and |6i| 7^ I62I, it is sufhcient to consider the 
case when —I<6i<0<62 = l. We call this case the “trapping magnetic 
step” (the term trapping will appear clearer after stating our results). 

( 3 ) If |6i| = I62I, we may restrict ourselves to the cases (61,62) = ( 1 , 1 ) and 
(61,62) = (— 1 , 1 ). This last one may be called the “symmetric trapping 
magnetic step” or “magnetic barrier”. 

( 4 ) If 61 and 62 have the same sign, we may assume that 0 < 61 < 62 = 1 . We 
call this case the “non-trapping magnetic step”. 

A physical description of the distinct classical orbits associated with these mag¬ 
netic field configurations is provided in [ 24 , 30 ]. The last two cases, ( 3 ) and ( 4 ), 
have already been studied in the mathematical literature (see the overview below). 
Our goal is to provide a mathematical treatment of the unstudied cases ( 1 ) and 
( 2 ). We will also give the consequences of our results on the behavior of quantum 
particles submitted to the magnetic field b. 


1.1. The band functions and their link with the magnetic Laplacian. 


• Direct fiber decomposition. In order to perform the spectral analysis, we can use 
the translation invariance in the y-direction and thus the direct integral decompo¬ 
sition (see [ 29 , XHI. 16 ]) associated with the Fourier transform with respect to y, 
denoted by J-y, 

TyLbT-^=\ ()b(fc)dfc, 

dfceK 


where 


t)b{k) = + Vb{x,k), with V\,{x,k) = {k — ab{x))'^. 


Note that the effective potential Vb{x, k) has the form 


Vb{x,k) 


{k — bixf, X < 0 , 
{k — b2xf, a: > 0 . 


( 1 . 2 ) 
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Figure 1 . The potential Vb{x,k) corresponding to case ( 1 ). 




Figure 2. The potential Vb{x,k) corresponding to case (2). 




Figure 3. The potential Vb{x,k) corresponding to case (3) with 
bi = —1 and 62 = 1. 




Figure 4. The potential Vb{x,k) corresponding to case (4). 


The domain of \)b(k) is given by 

Dom(l)(,(A:)) = {i/' G Dom(qt,(fc)) : {DI + Vbix, k))'ilj G L^(M)}, 
where the quadratic form qb(fc) is defined by 

[ \'tlj'{x)\'^+ Vb{x,k)\ij\'^dx. 

J 
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The spectrum of Cb, sp(£t,), satisfies 

spl-Cb) = U sp([}t,(A;)). 

feGffi 

We now define the n-th band function below the essential spectrum: 

Notation 1.1. We denote by Xb,n{k) the n-th Rayleigh quotient of We 

recall that if \b.n{k) is strictly less than the infimum of the essential spectrum, it 
coincides with the n-th eigenvalue of \)b{k). 

• State of the art. The cases (3) and (4) have already been extensively studied. 
Case (3) for which ( 61 , 62 ) = (Ijl) is just the Landau Hamiltonian. Its spectrum 
is pure point and consists of the so-called Landau levels; infinitely degenerated 
eigenvalues (2n— 1), for n > 1. The band functions are constants Xb,n{k) = 2n — 1. 

The case (3), for which ( 61 , 62 ) = (—1,1), can be linked to operators acting on 
a half-line. We first need to introduce the standard de Gennes operators, that will 
also be used in some results of this article. We recall below well-known properties 
of the spectrum of these operators. 

Notation 1.2. We let Eq = 0 and for n > 1, = 2n — 1. 

Notation 1.3. For o = N, (resp. o = D), let us introduce i)°{k) the Neumann 
(resp. Dirichlet) realization on K_|_ of D^+{t — k)^ and its eigenvalues {fig{k))i>i. 
For all £ > 1, the function admits a unique and non-degenerate minimum at 
k = denoted by Qe-i, with 0f_i € Moreover, ^^_i is also the 

unique solution of the equation ^^(/c) = k^. The function ^ decreasing 

function of k. Both ^J.f{k) and nf{k) go to -boo as fc —>■ —00 and go to E^ as 
k —>■ -boo. 

We refer mainly to [ 6 ] for the Dirichlet operator and [5, 13] for the Neumann 
one. All these properties can also be found in [10, 19]. 

Using symmetry arguments, it is proved in [9, 25] that for b = (—1,1) and i >1, 

Afa, 2 £-i(fc) = li^(fc) and Xb, 2 eik) = nf (k) keR, 
linking the de Gennes operator with the model with symmetric trapping magnetic 
steps. This analysis of the band functions provides a description of a quantum 
particle localized in energy far from the thresholds {En} and submitted to the sys¬ 
tem: they can exhibit edge current that may be viewed as quantum manifestations 
of the classical snake orbits [30] along the edge x = 0. These edge currents are 
exponentially decreasing far from the edge. The corresponding edge currents and 
the asymptotics of the band functions are described in [9]. The effective potential 
(see (1.2)) has the form of two symmetric double wells for fc > 0. The fact that 
each band function is exponentially close to a Landau level as k ^ -boo is called a 
splitting of the band function. The analysis of quantum states localized in energy 
near the thresholds is made in [16] for one edge quantum hall system and relies on 
a precise asymptotic analysis of the band function /i^(fc) as k —>■ -boo. The band 
functions associated to this case are displayed in Figure 5. 

The case (4), a non-trapping magnetic step, was treated in the initial paper by 
Iwatsuka [20] (along with the more general situation for which the magnetic held 
depends only on one direction and simply has constant values at ± 00 ). In this case, 
all the band functions are monotone nondecreasing. The analysis of the associated 
edge current is done in [17]. 
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Figure 5. Band functions for the symmetric double-well case b = (—1,1). 


Notice that all cases (l)-(4) are described from a physical point of view in [30]. 

1.2. Main results. We first consider case (1) for which ( 6 i,& 2 ) = (0,1). This was 
not treated in [18] or in [9]. It is physically interesting because the classical orbits 
of an electron originating in the half-plane a; < 0 are unbounded straight half-lines 
or lines (depending upon whether the orbit intersects a; = 0 or not), whereas those 
in the half-plane a: > 0 are bounded circles provided the Lamor orbit does not cross 
the edge at a; = 0 . 

Theorem 1.4. Let b = (0,1). 

If k < 0, then the fiber operator \)b{k) has only purely absolutely continuous 
spectrum [fc^,oo). Consequently, the spectrum of Lb is [0, oo). 

If k > 0, the essential spectrum o/f)t,(fc) is [k'^,oo). The band functions A(,,„ are 
nondecreasing functions with En-i < Xb,n{k) < En and lim/j_,.oo ^b,nik) = = 

2n — 1. For n > 1 and < k < (recall Notation 1.3), the operator t)[,(A:) 
admits exactly n simple eigenvalues below the threshold of its essential spectrum. 

The case (2) of a trapping magnetic step is technically more challenging. It 
relies on the asymptotics of the eigenfunctions of a Schrodinger operator with an 
asymmetric double-well potential. The results concerning the splitting of two energy 
levels depends upon the ratio of the two constant magnetic fields 61/62 = 61 since 
62 = 1 by normalization. 

Definition 1.5. We define £ = {En,n > l}[j{\bi\En,n > 1} and we define the 
splitting set 6 = ^ !}• We will say that there is a splitting 

of levels if 6 ^ 0 i.e. if there exist positive integers m < n such that 

2n — \ 

b=- - with6=16il. (1.3) 

Zm — 1 

As we shall see, the problem of the behavior for large k of the band functions 
can be reformulated into a semiclassical problem for a one dimensional Schrodinger 
operator with multiple wells (see for instance [14, 15]). In particular the limits 
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of the band functions are precisely the elements of £, counted with multiplicty. 
Roughly speaking, either one or two band functions converge toward each element 
of £. 

However our investigation will not rely on the semiclassical tools: We will use as¬ 
ymptotic expansions of special functions in order to derive the asymptotic behavior 
of the band functions. Let us now state our main theorem. 


Theorem 1.6. Let b = (5i,l), with —1 < bi < 0. For fc € M, the fiber opera¬ 
tor \]b{k) has only purely discrete spectrum Xb,n{k) with simple eigenvalues, and 
limfc_j._oo A(,,n(fc) = -boo. The set of limit points of the band functions as 

k —>■ -boo is precisely £. Moreover, we may describe the asymptotics of the band 
functions in the following way. 

(i) Non-splitting case. Let A S £ \ ©. We can write in a unigue way A = i?„ or 
A = bEn with n > 1. Then there exists a unique Pn € N* such that: 

(o-) A(,,p„(fc) — En = Enik) —)■ 0 (in the case A = En), 

k—¥-\-oo 

(b) Xb,p„{k) — bEn = En{k) —>■ 0 (in the case A = bEn). 

fc—V + OO 

In the first case, 

en{k) = -^2-2(1 + + o(l)), 

Vtt (’^-1)' 

and in the second case, 

e„(fc) =-+ b) , + o(l)). 

VTT (n-1)! 

(ii) Splitting case. Assume that the splitting set is not empty and consider A € ©. 
Let us consider n > 1 and m > 1 such that A = En = bEm. There exists 
Pn € N* such that 

j At,,p„(A:) En = £„ {k) 

1 Xb,Pn + ^{k) En = s'fi^ik') 


with, as k tends to -boo, 


£n (k) = - 


stik) = 


2"-3/2(i fc) 

(n - 

2™+l6-'"+l 

(to — 1)!(1 -b b)y/TT 




^2m+l 


'(1 + 0 ( 1 )), 

e-'=Vb(i + o(l)). 


The energy levels for which the exponentially small splitting occurs when A: —> oo 
behave like the band functions in the symmetric magnetic field case case bi = 
— 1 and 62 = 1. In that case (see [9]), two band functions approach the Landau 
level 2 n — 1 as fc —>■ 00, one is monotone decreasing and the other has a unique, 
nondegenerate minimum (see Figure 5). In the present article, with the help of our 
asymptotic expansions, we may deduce the existence of a minimum for some band 
functions. 


Corollary 1.7. With the same assumptions and notations as in Theorem 1.6, the 
following statements hold: 

(i) Let A € £ \ ©. Then each band function converging to A admits a global 
minimum. 

(ii) Let As©. Then k 1 —>■ Xp^(k) admits a global minimum. 
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Notice that the first band function always admits a global minimum. Moreover 
if |6i| is not a ratio of odd integers, all the band functions admit a global minimum. 
We conjecture that the only case when monotonicity occurs is Case (ii) of Theorem 
1.6 and corresponds to Ap^+i. In that case, the non-monotonicity of each band 
function depends on a non trivial way on the index of the eigenvalue. This situation 
seems striking to us, especially compared to other known case (see [9, 17, 34]). 

1.3. Numerical computations. We display on Figures 6-8 the band functions as¬ 
sociated with the magnetic step b = (—|, 1) together with the asymptotics provided 
below as fc —>■ -l-oo. The splitting occurs near each value in © = {2n — 1, n G N*}. 
These graphs have been obtained by using Mathematica. 



Figure 6. Band functions k ^ Xb.n{k) for b = (—|, 1). 
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Figure 7. Zoom near A = | on the band function k i-A \, 4 {k) 
for b = (—^,1). In dash line: the asymptotics from Therorem 1.6, 
case (ib). Here n = 3. 



Figure 8. Zoom near A = 1 on the band functions k Xb,p{k) 
for b = (—1,1) and p = 2,2, (splitting case). In dash line: the 
asymptotics from Therorem 1.6, case (ii). Here A = 1, n = 1 and 
m = 2. 


1.4. Consequence of the analysis of the band functions. In this section we 
qualitatively present possible consequences of our results. 

• Standard perturbation of the Hamiltonian. If we perturb the Hamiltonian by 
a suitable electric potential, discrete spectrum might appear under the essential 
spectrum, which is the infimum of the first band function. In particular, the exis¬ 
tence of a minimum is useful if one wants to study the counting function of created 
eigenvalues. Here we know that the infimum is reached for some finite frequency 
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k*. Then, following the standard procedure described for example in [27], the num¬ 
ber of discrete eigenvalues can be described, depending on the degeneracy of the 
minimum and on the decay of the electric perturbation. 

In the same way, it is possible to derive a limiting absorption principle at energies 
close to the bottom of the spectrum of Cb (and also near each band minimum), 
as in [33]. The case of thresholds corresponding to the limits of band functions 
can also be studied in the spirit of [26]. In particular the asymptotics of the band 
functions for large k allow to describe the suitable functional spaces when writing 
a limiting absorption principle near the Landau levels. 

• Semiclassical magnetic Laplacian with piecewise constant magnetic field. Let us 
consider the semiclassical magnetic Laplacian {—ihV -I- A)^ in the plane with a 
magnetic field such that 


V X A = B = 6iln -l- 62lcnj 

where 12 is a smooth and bounded domain. It is well-known that the semiclassical 
analysis of the low lying spectrum (motivated by the analysis of the third critical 
field in the Ginzburg-Landau theory) is governed by the spectral behavior of some 
local models (see [10, 3]). When the magnetic field is not continuous, models 
operators in the form of £(, have to be considered. It is known that the minima 
(which are often non degenerate) of the corresponding band functions may combine 
with the curvature (of the boundary for instance) to generate the semiclassical 
asymptotics. In particular, the existence of a minimum of the band functions 
implies a microlocalization of the magnetic eigenfunctions (see [10, 28, 4]). 

• Bulk states localized in energy near the thresholds. In general Iwatsuka Hamil¬ 
tonians {i.e. Hamiltonians associated with magnetic fields depending only on one 
variable), the analysis of quantum states submitted to the magnetic field and lo¬ 
calized in energy in an interval I depends on the behavior of the band functions 
when crossing the energy interval I, see [21]. For example if / does not contain 
any threshold, then the quantum states are localized in space near the variations of 
the magnetic field (here, the jump line {x = 0}) and their current (defined below) 
is not small. This corresponds to the fact that a classical charged particle moves 
along the edge {x = 0}, in a “snake orbit”. 

One may also show, as in [6, 16], that there exist states, called bulk states, 
localized in energy near the thresholds, bearing small current. If the threshold 
corresponds to a limit of a band function, these bulk states are exponentially small 
near the jump line. The precise (de)localization of these states and the current they 
bear can be derived from the asymptotics of the band function by the following 
method (see [16]). Assume that I = {X — 6, X + 6) where A G 6 and <5 | 0. 
The current of a state if G L^(M) along the jump line is given by (Jyipjip) where 
Jy := ^[Cb,y] = —idy — ab{x) is the current operator. This quantity may be linked 
to the derivative of the band functions by using the Feynman-Hellmann formula. 
Using the asymptotics from Theorem 1.6, one may show that the derivative of the 
band functions is small on the large k interval corresponding to the values where 
Xn{k) G I. This would provide an upper bound on the current carried by bulk 
states, depending on how far from the threshold is their energy. Then follows an 
analysis in phase space of these states in order to show that delocalization in the 
frequency k involves delocalization in the a;-variable. 
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Note finally that the sign of the derivative of the band function is related to the 
direction of propagation of the particle along the jump line of the magnetic field. In 
particular quantum states related to non-monotonic band functions (see Corollary 
1.7) may flow along the line in two possible directions. This corresponds to the 
possible negative velocity described in [30]. We refer to [34] for an analysis of the 
direction of propagation of quantum states in a translation invariant magnetic field. 

2. Magnetic wall 

This section is devoted to the case b = (0,1). Away from the edge at x = 0, a 
classical electron moves in a straight line for a: < 0 and in a closed circular orbit 
for a: > 0. If the orbit intersects the edge, then the classical electron escapes to 
X = —oo. The absence of magnetic field for x < 0 results in a half-line of essential 
spectrum for the quantum electron. As we will see, the spectrum varies with the 
sign of k. 

The reduced potential 14 satisfies I4(x, fc) = for x < 0, and t4(x,fc) —>■ -boo 
as X —>■ -boo. As a consequence, we easily deduce the following proposition. 

Proposition 2.1. Assume that b = (0,1) and k G M.. The essential spectrum of 
\]b{k) is given by 

sPe5s(f)6(fc)) = [A:^-boo). 

Moreover, sp(£(,) = [0,-boo). 

In fact, we can prove slightly more. 

Proposition 2.2. Assume that b = (0,1) and fc S K. The operator i)b{k) has no 
embedded eigenvalues in its essential spectrum. 

Proof. Let us consider A > and ip G Dom(()(,(fc)) such that: 

- Ip” + {k - ab{x)yip = \ip. (2.4) 

For X < 0, we have —ip” = (A — k‘^)ip whose only solution in L^(M_) is zero. 
But since the solutions of (2.4) belongs to and are in C^(IR), this implies that 
■0(0) = ip'{0) = 0 and thus -0 = 0 by the Cauchy-Lipschitz theorem. □ 

Let us now describe the discrete spectrum, that is the eigenvalues A < Since, 
for fc < 0, we have qb(fc) ^ we deduce the following proposition by the min-max 
principle. 

Proposition 2.3. Assume that b = (0,1) and k < 0. Then 
spii)b{k)) = spg, 3 (f)[,(fc)) = [fc^-boo). 

Therefore we must only analyze the case when fc > 0. The following lemma is a 
reformulation of the eigenvalue problem. 

Lemma 2.4. Assume that b = (0,1) and k > 0. The number X < k'^ is an 
eigenvalue of f)b(A:) if and only if there exists a non zero function ip G L^(K_|_) 
satisfying 

j —ip”(x) -b (x — k)'^ip = Xip{x), (2.5a) 

\ ip'{0) - Vk^-XipiO) = 0. (2.5b) 

Moreover, the eigenfunctions o/f)b(fc) can only vanish on K+. The eigenvalues of 
\)b{k) are simple. 
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Proof. We consider = Xip. Since ip G L^(R-), there exists A G K such that, 

for a; < 0, ip(x) = Then we have to solve —ip" + Vb{x,k)ip = Xip for 

X > 0 with the transmission conditions 

i/>(0) = A and ip'(0) = A\J— X 

or equivalently 

ip\0)-^/k^^'ip{0)=0. 

In particular, A cannot be zero. The simplicity is a consequence of the Cauchy- 
Lipschitz theorem. □ 

Lemma 2.5. Assume that b = (0,1). The functions M+ ^ k i-G are non 

decreasing. 

Proof. We use the translation x = y -\- k to see that f)6(fc) is unitarily equivalent to 
+ V{y,k) with V{y,k) = t(-oo,-k){y)k'^ + l(-fc.+oo)(2/)y^- For 0 < fci < fe, we 
have: 

V{y,k 2 ) -V{y,ki) = + l(_oo,-fc 2 )(y)fc| - l(-oo,-fei)(j/)fc? 

— ^(—fc2,—fei)(y)^i "f i(—oo,—fc2)(2^)^2 i(—oo,—fci)(y)^i 

= l(_oo,-fc2)(y)^2 - l(-00,-fe2)(2/)^l > 0- 

By the min-max principle, we infer the desired monotonicity. □ 

The next lemma is a consequence of the Sturm-Liouville theory (see for instance 
[32]). 

Lemma 2.6. Assume that b = (0,1), and let n G N be such that A[,^„(fc) < 
Then, the corresponding eigenspace is one dimensional and is generated by a nor¬ 
malized function ip\,,n{k), depending analytically on k. Moreover, ipb,n{k) has ex¬ 
actly n — 1 zeros, all located on the positive half line. 

Proof. We have only to explain the part of the statement concerned with the zeros. 
Thanks to Lemma 2.4, one knows that the zeros are necessarily positive. Let us 
briefly recall the main lines of the Sturm’s oscillation theorem. Let us denote by 
m — 1 the number of zeros ot tp = ipb,n{k) and let us check that m < n. For 
j = 1,... ,m — 1, Zn,j denotes the j-ih zero of ipb,n(k) (the zeros are simple and 
isolated by the Cauchy-Lipschitz theorem) and we let pj = -^Zr, j+i) extended 

by zero. By definition, 

= / ip'jLp'iPVb{x,k)ipjg}(,Ax, 

JM. 

SO, 

JM. 

We apply the min-max principle to the space Fm = span ipj to get Xb^m < ^b,n 

and thus m < n. In particular we recover that ipb.i does not vanish. Let us explain 
why ipb ,2 vanishes exactly once. Introducing the Wronskian ITi _2 = 'f’i'4’2 ~ V'2'0ij 
we notice that IT[ 2 = (Ab,i — At,, 2 )i/’b.iV’b, 2 - If V'b .2 > 0, then IT[ 2 < 0- Since the 
eigenfunctions belong to the Schwartz class at infinity, Wi^ 2 {x) —>■ 0 as |x| —)■ -boo. 
This is a contradiction. With the same kind of monotonicity arguments, we can 
prove that the zeros of the eigenfunctions are interlaced and that ipb,n has exactly 
n — 1 zeros. □ 
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By using the harmonic approximation in the semiclassical limit (see [7] and the 
Section 3.1 below), we can prove the following lemma. 

Lemma 2.7. Assume that b = (0,1). For all n> 1, we have 

lim = En. 

k—¥-\-oo 

In particular, for k large enough, we have Ab,n(^) < En < . 

Let us now prove that the n-th band function lies between the two consecutive 
Landau levels En-i and E^. 

Proposition 2.8. Assume that b = (0,1). For all n > 1 and for all k > 0 such 
that Ab,„(/c) < we have Xb.n{k) € 

Proof. By Lemmas 2.5 and 2.7, A(,,n(fc) < (the strict inequality comes from the 
analyticity). It remains to prove that Xb,n{k) > En-i- Clearly Ai,^i(A:) > Eq. We 
rely on Notation 1.3. For n > 2, we consider the function ipn{x) = f’b,n{A + Znp{k)) 
which satisfies l)^(fc —= At,_„(fc)(/?„ and has exactly n—2 zeros on K+, cf. 
Lemma 2.6. By the Sturm’s oscillation theorem, is the {n — l)-th eigenfunction 
of [}°(fc — Zn,i{k)). Therefore we have Xb.n{k) = fi^-iik — Zn,i{k)). Moreover for 
all € > 1 and fc G K, g^{k) > Eg. This follows by domain-monotonicity, once we 
notice that [}'^(fc), after a translation, is given by the harmonic oscillator 
on {—k, -foo) with Dirichlet boundary condition at a; = —k, together with the fact 
that the harmonic oscillator D^+x"^ on K has simple eigenvalues Eg. This provides 
the desired conclusion. □ 

We recall that the sequence (^n)n>i is defined in Notation 1.3. For the next 
proposition, we will need to use the so-called parabolic cylinder functions. We 
denote by U (a, x) the first Weber parabolic cylinder function (whose properties 
are summarized in the Appendix, section A; for more details see, for example, [1, 
section 12] or [2]) which is a solution of the linear ordinary differential equation: 

- y"{x) + ^x'^y{x) = -ay{x). (2.6) 

The Weber function U{a,x) decays exponentially when x —)■ +oo. 

Proposition 2.9. Assume that b = (0,1). For all n> 1, the eguation Xb^n{k) = k^ 
has a unique non negative solution, k = ^n-i; such that, locally, for k > kn, 
^b,n{k) < Moreover we have Ab,n(Cn-i) = ©n-i- 

Proof. Thanks to Lemma 2.7, we can define kn = max{/c > 0 : Ati,n(fc) = fc^}. By 
continuity, we have, for all k > kn, A(,,n(fc) < k'^. Let us now prove the uniqueness. 
Let us consider a solution kn > 0. For all integer p > po with po large enough, we 

have A(,,n (j^n + p) < {kn + . Let us now consider the eigenvalue equation 

T ikn,p ^b(x)} Fn,p ~ ^n,pFn,p, 

where ipn,p = tpb.nikn.p), fcn,p = kn + ^ and Xn,p = Xb,n{kn,p)- Let us investigate 
the limit p -l-cx). As seen in the proof of Lemma 2.4, we know that there exists 
a G K* such that, for a: < 0, 


/ \ x\[y^ 

= ae V - 


-A, 
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By solving (2.7) on a; > 0 and using the parabolic cylinder function U, we find that 
there exits /3 such that, for x > 0, 

^n,p{x) = PU - kn,p)^ ■ 

Since ^ 0, we have {a, P) ^ (0,0) and (pn^p is at x = 0, we get the 
transmission condition: 

\/~klp-KpU -~kn,pVi^ - V2U' -~kn,pV2^ = 0 . 

By continuity and taking the limit p —>■ +oo, we get 

U' =0. 

Notice that the function x ^ U •\/2(x — fc„)^ solves the differential equation 

I - y"{x) + (x - knfy{x) = kly{x) 

\ 2 /'( 0 ) = 0 and ?/(0) ^ 0 

Moreover it belongs to 5(M+). Therefore, there exists ^ > 1 such that pf {kn) = 
and therefore kn = and = 0^_i. By Proposition 2.8, we know that 
kn = ^b,n{kn) G [En- 1 , En]. Moreover, we recall that 0f_i G {Ei-i,Et). This 
implies that i = n. □ 


3. Trapping magnetic step 


In this section, we focus on the case b = (6i, 1) with —1 < 6i < 0. The case of 
6i = — 1 is treated in [9]. In that case, due to the symmetry, all levels are split as 
k —> + 00 . The asymptotic expansion comes from [16, 25]. The band functions for 
this case are shown in Figure 5. 

When bi < 0, the effective potential Vb{x;k), described in (1.2), diverges to 
infinity as x — >■ ±oo, so that for all fc G K, the operator f)b(fc) has compact resolvent. 
The spectrum of f)b(fc) is discrete and the eigenvalues (At,_„(A:))„>i are all simple. 

For fc < 0, the effective potential 14 (x; k) has the form of an asymmetric potential 
well, see Figure 2. It follows that Vt,{-]k) > so that, for fc < 0, A(,,„(fc) > k'^. In 
particular, we have, for all n > 1, 

lim Ati,n(fc) = + 00 , k < 0. 

fc—>• —OO ’ 


For fc > 0, the effective potential has the form of an asymmetric double-well: 


Vb{x;k) 


{k — 6ix)^ X < 0 
{k — x)^ X > 0 


(3.8) 


This asymmetric double-well admits two nondegenerate minima at x = ^ and 
X = fc, for fc > 0 and hi < 0. 

Let us now investigate the behavior of the band functions in the limit k 
+00 and prove Theorem 1.6. Due to the asymmetry of the double-well potential 
(3.8), we will prove, among other results, that not all the Landau levels exhibit an 
exponentially small splitting as in the case of {hi, 62 ) = (~ 1 ) !)• 
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3.1. Limits of the band functions. Let us first describe the set of the limits of 
the band functions At,,„(A:) as fc —>■ +oo. As noted above, the limit as fc —> —oo of 
any band function is plus infinity. 

Lemma 3.1. Assume that b = (bi, 1), with —1 < bi < 0. Recall from Definition 
2.1 that £ denotes the union of the Landau levels for magnetic fields bi € (—1,0) 
and 62 = 1; O'lT'd that © denotes the splitting set {Ej^^n > l}n{|5i|i?„,n > l}.Then 
the set of limits as k ^ +00 of the band functions is precisely the set of Landau 
levels C: 

( lim Afa,„(fc),n > 1 } = £. 

Ifc—V+OO J 

The band functions exhibit the following behavior as k ^ +cx); 

(1) Non-splitting case. If a band function limit A € £ \ 6 , then for any 

eg > 0, there exists a K > 0 such that if k > K then it holds that 
rankl[>,_£j, = 1 and there exists a unique n > 1 such that 

A[),n converges towards A. 

(2) Splitting case. If a band function limit A S ©, then for any eg > 0, there ex¬ 
ists aK > 0 such that ifk > K then it holds that rank = 
2 and there exists a unique n > 1 such that Xb,n and A(,,n+i converge to¬ 
wards A. 

Proof. Let us first provide a semiclassical reformulation and summarize the stan¬ 
dard procedure of the harmonic approximation for multiple wells (see [12, Chapter 
4] for example). For k > 0, we let x = ky and the operator l)b(fc) is unitarily equiva¬ 
lent to where l)b,h is the self-adjoint realization on L^(M) of the differential 

operator h'^Dy -\- Vb{y; 1). The minima of Vt,(s 1) are reached for x = bf^ < 0 and 
X = 1. We have Vf,"(l; 1) = 2 and Vf,"( 6 ["^; 1) = 2|&i|. Let us just recall the spirit 
of the analysis. Let us consider 77 > 0 such that bf^ -I- ?7 < 0 and 1 — ry > 0. Then, 
let us consider the operator ^b,h,r] = Tb,n,ri ® where and are the 

Dirichlet realizations of hfDy + V(y, 1) on L^((— 00, 61 -I- 77 )) and on L^((l — 77 ,-l-oo)) 
respectively. Let us now consider an arbitrary constant C > 0. It is well-known 
(see Helffer’s lecture notes [12]) that any eigenfunction of ^b,h,Tj or of {)b,h asso¬ 
ciated with an eigenvalue X < CH satisfies decay estimates of Agmon type away 
from the minima. In particular, this allows to prove the existence of c > 0 such 
that sp(^t,,s) n {A S K : A < Ch} and sp(^b,?i, 77 ) n {A S M : A < Ch} coincides 
modulo 0[e~'^^ ). Then, by the harmonic approximation, the spectrum of i)b,h,r) 

is well-known in the semiclassical limit: the 7r-th eigenvalue of '^bhr] equivalent 
to \bx\Enh and the one of i)b.h,rj ^o Enh. The desired conclusion follows. The last 
part of the statement follows from the fact that the spectra (with multiplicity) of 
f)b,h and i}b,h,Ti coincide modulo an exponentially small error. □ 

As described in [16], the convergence of the band functions toward these limits 
corresponds to the existence of bulk states with energies close to the Landau levels. 
We next provide a more accurate description of the convergence of the band func¬ 
tions towards their limits. This is important for two reasons. First, this allows a 
more quantitative characterization of the associated bulk states following the anal¬ 
ysis of [16]. Second, the manner in which the band function approaches its limit 
as fc —^ -koo, that is, whether it approaches from above or below the Landau level, 
provides a hint as to whether the band function is globally monotonic or not. In 
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order to describe the asymptotic behavior of the band function, we will leave the 
semiclassical techniques and use the characterization of the eigenvalues Af,,„(fc) as 
zeros of certain special functions. 

3.2. Splitting and non-splitting of levels. We consider k > 0 and calculate 
the behavior of the energy levels in the splitting and the non-splitting cases. This 
is a refinement of Lemma 3.1, where we have shown that each band function 
tends to a threshold that is an element of the set of Landau levels £. There is no 
explicit formula for the index of the corresponding Landau level as a function of n, 
the band index. 

In this section, for the sake of simplicity, we will consider a converse strategy: 
Given a threshold or bEn, where n G N* is fixed, we describe the asymptotic 
behavior of the band functions converging toward this threshold. As described in 
Lemma 3.1, there are either one or two band functions Xb,pni^) converging toward 
a threshold En or bEn in £. Although the indices of these band functions are 
not explicit in terms of the index n, the asymptotics of the related band functions 
A[i,p„(fc), as k ^ -bcxD, are explicit functions of n. 

Let us notice that this procedure will provide the asymptotics of all the band 
functions since there is a bijection between their limits and the set of the thresholds 
£ (counted with multiplicity). 

3.2.1. Algebraic characterization of the eigenvalues of \)b(k). We begin with an 
algebraic characterization of the eigenvalues of ^b{k). Recall that U{a,x) is the 
first Weber parabolic cylindrical function introduced in Section 2 as a solution of 
( 2 . 6 ). 

Proposition 3.2. Assume that b = (6i,l), with — 1 < 6i < 0. A number A is an 
eigenvalue of f)ti(fc) for k gM. if and only if it is a solution o/I?i,(A, k) = 0, where 

DbiX, k) := C/(- A, -V2k) 

+ VbU'{-E^,-V2f^)U{-l,-V2k) 

and b = |6i|. 

Proof. From the eigenvalue equation and the form of the potential (3.8), we are led 
to solve for X > 0: 

-y''{,x) + {x- kfy{x) = \y{x) 

with the condition that the solution is in L^(M+). The only such solution is 

X !->■ [/(—|, •\/2(x — k)). 

In particular, when A = 2n — 1 is a Landau level, the solution is the n-th Hermite 
function centered at x = k. On the half line x < 0, we solve 

-y"(x) -f {bix - kYy{x) = Xy{x) 

with the condition that the solution is in L^(M_). This solution may be expressed 
as 

x^U{-^,V2{-Vbx- ^)). 
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Therefore, if is an eigenfunction associated with A = Ab,n(fc)) it may be 

written as (up to an overall multiplicative constant): 

"'-Aif-o 

{ aU{—^,V2{x — k)) if cc > 0 

where a S K is a constant. Since the eigenfunction is at 0, we have 

- = aV2U'(-l,-V2k) 

It is possible to solve the above conditions if and only if the 2x2 matrix 

/ [/(-A,_y2^) u{-l,-V2k)\ 

u'{-^,-V2k)) 

is non-invertible. □ 

We will use the following result on the asymptotic expansion of the determinant 
Db{X, fc) as fc —)• + 00 . 

Lemma 3.3. Assume that b = with —1 < 5i < 0, and let b = |6i|. The 

determinant Db{X,k) of Proposition 3.2 can be written as 

Db(X,k) = e(i+E)Tri(A,fc) + eAi+i)TT2(A,fc) 

+ e(i-5)^T3(A,fc) + e(-i-^)^r4(A,fc), 

where the factors Tj(X, k) admit the following asymptotic expansions as k ^ +oo.' 

Ti{X,k) = ^^^ 

21 
1 


(3.9) 


1 b + 1 \ 3 {l + b)A b+A / ^ n, \ 

3^fc“^^22 ib— TT (l + 0{k~^)) , 


r(^)r(V) 


A 2 (A, k) = (1 + b){- sin ^) (1 + 0{k-^)) , 

1 


nix,k) = 


r(V) 


-k 26 


^A-2 


2^ 46 " 6 46 (1-I-5)(_ sin ^) (l + 0(fc ^)) , 


TiiX, k) = V2k-^^2^'^b-^ cos (1 + 0(fc”^)) ■ 


2b 

Proof. These expansions follow by the asymptotic expansion of the Weber function 
U. We refer the appendix, section A, and to [2, Chapter 19]. □ 

Remark 3.4. We have controlled these rather tedious calculations with the computer 
software Mathematica. 

Remark 3.5. To obtain the large k asymptotics of Db(X, k), we factor out the leading 

exponential factor from the right side of (3.9). The zeros of Db{X, k) are 

the same as the zeros of 


Db{X,k) := e-(i+^)^I4a(A,fc), 


where 


Db{X, k) = ri(A, k) + e-'=V2(A, k) + e-TT3(A, k) + e-(i+5)'=V4(A, k). (3.10) 
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We next note that the functions Ti and T 2 are bounded as A: —>■ + 00 . The coefficient 
T 3 grows like for A large enough. Similarly, the term T 4 is increasing in 

k. However, the prefactors of these terms exhibit Gaussian decay. 

The strategy of the proof of Theorem 1.6 consists of writing the band function 
Xb,p„{k) as 

where limfc ^+00 A[,,p„ (fc) = En, with either En = En or En = bEn, and limfe_j.+oo Sn{k) = 
0. For all fc G M, we have 

Db{Xb,p^{k),k) = Db{En +£n{k),k) = 0. 


The distinction between the non-splitting and splitting cases is encoded in the 
gamma function coefficients of the Tj(A, k) terms in Lemma 3.3. 

In the non-splitting case, we have for one of the gamma functions 




^ T 



In case (ia), the argument of the gamma function is ((6-1- l)/2 — n)/b that is never 
a negative integer so T”^ does not have a zero and does not contribute to the 
asymptotics. In case (ib), the argument of this gamma functions is (1 — n) so there 
are poles of the gamma function. The corresponding zeros of T”^ contribute to the 
vanishing rate of Ti and T 2 . As for the other gamma function, we have 


1 Xb,pni.k) 




r (1 — n) (ia) 
T{^-nb) (ib) 


so there is a pole in case (ia). These asymptotics are developed in Section 3.2.2. 

In the splitting case, we have that b — En/E^, and it follows that both gamma 
functions have poles. Consequently, both factors T”^ have a zero as A: —>■ -|-oo and 
these contribute to the asymptotics as developed in Section 3.2.3. This analysis also 
shows that Ti{En,k) = 0 since for (ia) the inverse of the gamma function r(i^) 
vanishes for A = E^ and in the case (ib) the inverse of the other gamma function 
r(^)=0for A = 6 F;„. 


3.2.2. No splitting of levels. Let us prove Item (i) in Theorem 1.6. The existence 
of Pn comes from Lemma 3.1. Let us assume case (ia), the proof of case (ib) being 
symmetric. Then 


lim 

k —^ “|“Oo 


r( 


b-A,, 


dfc) 


26 


) r(^) 


^ 0 , 


(3.11) 


the last relation coming from the facts that the only roots of l/T are the negative 
integers and that ^ —N, since this is not a splitting case. Recall that the 

roots of l/T are simple. In particular, as k gets large. 


r( 


l-Ai, 




( 1 )' 0(e„(A:)) 

1 ) 1(1 + 


(3.12) 


As mentioned above, we replace At,,p„ (k) by A„-|-e„(A;) in the equation DbiXb,p„ (k), k) = 
0 and obtain an equation in the form 

a{k)£n{k) + P{k) = 0 , 














18 


P. D. HISLOP, N. POPOFF, N. RAYMOND, AND M. P. SUNDQVIST 


where a and /3 are analytic functions such that 


a{k) 


k—>-\-oo 




' ■n/' b—En \ 
^ 1 2& 1 


i>+l p 3 

-k 26 -®" 22 ' 


(l + t)Er. 


h + Ejt, 

b “•>> IT, 


m 


k—^+OD 


r(^) 


-k St 


-E„-2 


-)E, 


6-1 b+En , 

46 & 46 (1-1- sin 


Er,n 


,-fc 


obtained by using Lemma 3.3. We deduce the asymptotic behavior of en(fc) as 
k —>■ +00. 

Since e„(fc) is negative the corresponding band functions approach the 

Landau level En or bEn as fc —>■ +oo from below. Taking into account the fact 
that Xb,p^(k) —>■ +00 as k ^ —oo, this means that these band functions are not 
monotonic and have at least one minimum. 


3.2.3. Splitting of levels. Let us finally give the proof of Item (ii) in Theorem 1.6. 
When the condition for splitting is achieved, two band functions have the same 
limit as fc —> +oo. We prove that one level approaches the common Landau level 
from below, and hence has a minimum, whereas the other band function tends to 
its limit by above. 

The existence of is a consequence of Lemma 3.1. We also know from this 
lemma that both functions tends to 0 as fc —>■ +oo. We use the generic notation 
A(fc) = En + e(fc) satified by these two functions. Recalling that b = we note 
that as fc tends to +oo: 


and 



■^(f) 


Substituting these Taylor expansions into (3.9), the condition Db{X{k),k) = 0 
implies that are the roots of the equation 


Ae^+s{B + C) + D = 0, 


(3.13) 


where {A, B, C, D) are analytic functions of fc respectively equivalent to the follow¬ 
ing functions as fc —?► -foo: 


^ I 1 77T. T b-|-l 

A{k) = (-l)"+'”(n - l)!(m - l)!7r2-5- 

B{k) = (-1)"+™(to - 1)!v^ 2-3/2+^^^^^6^-^(1 -f b)k'^^^-‘^e-^\ 

C{k) = (-l)"+™(n - 1 )!v^ 2-3/2+^^^^S^6^-1(1 fc)fc^^"-2e-'='/'', 

D{k) = (-1)"+— 

(3.14) 

The discriminant 

A = (R -h (7)2 - AAD 
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is positive for k large enough since there hold B{k) 3> C{k) and B{k)‘^ D{k) as 
k goes to +00 (notice also that A{k) = o(l)). In particular, the two roots of (3.13) 
admit the following expansions: 


D 


B 


= —^(1 + 0 ( 1 )) and e =—^(l + o(l)). 

B A 

Therefore, using the equivalents given in (3.14), we get as k gets large: 

^-m+§2m+f 


£+(/c) = 


^2m+lg-fe 


(m - 1)!(1 + b)y/n 

e-{k) = ^ 

By definition, e~ < s^. Due to the obvious sign of e+ and e“ for k large enough, 
we get and the result follows from the above asymptotic expansions. 


Appendix A. Parabolic cylinder functions 


We list some properties of parabolic cylinder functions, and refer to [2, Chap¬ 
ter 19]. Given a e K, denote by U{a,x) and y(a, a;) the so-called standard solutions 
of the differential equation 

-y'ix) + Ix'^yix) = -ay{x), a: G K. 

If follows by the series expansions of U and V that they are real-valued if o G M 
and a; G K. Moreover, as a: —>■ -l-oo, 

(- -I- a) 

U{a,x) ^ exp(-crV4)cr-°-V^^(-l)^ ^^ (A.15) 

j=o ^ ’ 

and 

VM -^exp(a:V4)a:°-V^X] ^/!(2^2)y - (^-16) 

Here {x)j = x{x — l)(a; — 2) • • • (a; — j + 1) denotes the Pochhammer symbol. It also 
holds that 

ttV (a, x) = r(a + 5){sin(7ra)C/(a, x) + U (a, —a:)}. (A.17) 

This enables us to give an asymptotic expansion of U{a,x) as a; —>■ — 00 . Indeed, as 
X -boo. 


U (a, —x) 


r(a -b 4) 


+00 (1 


exp(A/4)x--/»^|^ 

j=0 ’ 


— sin(7ra) exp(—a;^/4)a; 


+00 (1 a) 

1=0 


2j 


j!(2a;2)a 


If we keep only two terms in the sum, 

U{a,-x) - exp(a;V4)a;°~^/^(l-b 4(4 - 0)24 +^(1A'^)) 

r(a -b |) \ ^ ^ x^ / 

- sin(7ra) exp(-a;^/4)a;“““^/^^1 - 4(a -b |)2^ + C>(l/a:^).^ . 


(A.18) 


(A.19) 
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Next, it holds that 


and thus 


d X 

U'(a,x) := —U(a,x) = —U(a,x) — U(a— l,x), 
ax 2 


U\a, —x) = — — U{a, —x) — U{a — 1, —x). 


Combining this with (A.18), we have, as x —>■ +oo. 


U'{a,-x) - r/ 1 ' exp(xV4)x°+^/^x 

r(a+ 2) 


1 (- 

1 V2 


2 ^ j!(2x2)J 

j=o ^ ^ 




^—' j! (2x^)1 

j=o ^ > 


+ sin( 7 ra) exp(—x^/4)x “+i/2x 


(«+5)2, 


’ j!(2x2)J ’ 


2j 


1=0 


1=0 


j!(2x2)J 


With two terms in the sum only. 


U'{a, —x) ~ — 


r(a + 2) 


exp(x^/4)x“''‘^/^x 


(2 { 4 ^^ ~ + (a - \)]^ 

+ sin(7ra) exp(—x^/4)x““''"^/^x 
(“2 + {-;^(a+ 5)2 + 5 ) 2 }^ + 0 (l/x'‘)). 


(A.20) 


(A.21) 
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